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1 Introduction 

In a recent paper [ 1] a complete list of exactly solvable cases was presented for 
the three-state vertex model with ice rule and certain symmetries. As a criterion for 
integrability the "additive" Yang-Baxter equation was used and solved exhaustively 
for the Boltzmann weights in the considered class of models. The list in [1] comprises 
several well-known models, but also four non-trivial new ones. The solution to these 
models was not given in [1]. It is the purpose of the present paper to present the exact 
solution of the four models, namely the calculation of the thermodynamic properties, 
i.e. partition functions, correlation lengths, critical exponents, etc. At the same time 
we solve the associated quantum spin-1 chains. 

* Work performed within the research program of the Sonderforschungsbereich 341, Koln-Aachen- 
Jiilich 

^Permanent address: L.D. Landau Institute for Theoretical Physics, Kosygina str. 2, GSP-1, 
117940, Moscow, Russia; Email: matveen@cpd.landau.free.net 
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The vertex models considered in this paper are defined on a square lattice where 
spin variables are placed on the bonds. Each spin may take three values, say or 
±1, and there are interaction energies associated with each vertex depending on the 
local spin configuration. The corresponding local Boltzmann weight is denoted by 
for a spin configuration (j,, u, a, and (5 on the left, right, lower, and upper bond 
of the vertex, respectively. 

The partition function for a lattice of size NxL with periodic boundary conditions 
is given by 

Z = TrT L , (1) 
where the transfer-matrix T is the product of the Boltzmann weights in a row: 

rpai,...a N \ ^ pf l "l pf 2 "2 pfHV «JV /n\ 

1 /3 1 ,...f3 N — 2^ /3i @2 ' ' ' m Pn ' \ Z ) 

fii,.../j, N 

As in [1] we impose the ice-rule 

a + fi = [3 + v (3) 

which leads to nineteen allowed, i.e. non-zero, Boltzmann weights, while the further 
symmetries 

f?M Q _ p-ju-a _ pct^t _ jdu/3 ( a\ 

n vf3 — n -v-(3 — n (3v — n fj,a- K^J 

reduce them to only 7 independent weights: 

a, = b = c = 

„ _ pio _ pl-1 „ _ pOl 

d = RZ (5) 

In the following we shall consider the four new models which in [1] have been 
derived as the solutions #8, #9, #2, and #3 of the Yang-Baxter equation and which 
we label as I to IV. The four models are defined by their respective (not normalized) 
Boltzmann weights which all depend on the spectral parameter u. Models I and 
II have no further variable interaction parameter and are parametrized (with slight 
changes as compared to [1]) as follows. With U = expu we have the weights: 

Model I: 

4 - U 4 2(U 4 - 1) _ U 4 + 2 

~W 2 ~ 1 ~ 3U~ 2 ' C ~ 3U 2 ' 

e = Tj' 9 = ~^2 b ' p = ' d = 1 ' ^ 



a = 



Model II: 



a(a 2 -U 4 ) , TT2 U 4 -1 a 2 U 2 
a = — -r , b = all 2 



U 2 {U 4 + aY U 4 + a U 4 + a 
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e = Ua, g = —b, p = 0, d = 1, (7) 
where ct = v/ ^ +1 . 

In the physical region all Boltzmann weights must be positive which leads to 

< Rew < A, A = {| nv l' f rm °f!!; (8) 

[ In y/a, for model II. 

The models III and IV depend on an additional interaction parameter besides the 
spectral parameter u and are parametrized differently in different regions. 

Model III: We have b = g = and a = c = d. With 

A = a2 + P 2 ~ e \ (-00 < A < oo) (9) 

the parametrization depends on the value of A. 
• A > 1: 

sinh(A + u) sinhw 

a— — , p— — — , e = 1, A = cosh A, (10a) 

smh A smh A 



• A < -1: 



sinh(A — u) sinhw A , , , „, , 

a = ^— — -, p= — — , e=l, A = -coshA, (10b) 

smh A smh A 

• -1 < A < 1: 

sin(A — u) sinw 

a = , p— -, e = 1, A = — cos A, (10c) 

sin A sin A 

with A > 0. 

Model IV: Again we have b = g = and a = c. However, now we require 

ad + p 2 - e 2 = 0, (11) 

which will prove to be a free fermion condition. Otherwise we shall use two different 
parametrizations 

• a) 

a = cosh u ± cosh A sinh u, d = cosh u =F cosh A sinh u, 

p = sinhAsinhw, e = 1, (12a) 

• b) 

a = cos u ± sinh A sin u, d = cos u =p sinh A sin u, 

p = coshAsinw, e = 1, (12b) 
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where A > 0. 

This completes the parametrization of the four models. We emphasize that all 
parametrizations satisfy the "standard initial condition" 

R^(u = 0) = 5^5 va . (13) 

One of the consequences of the Yang-Baxter equation is the existence of a family 
of commuting transfer matrices T(u) generated by the spectral parameter u 

[T(u),T(v)]=0. 

As a result of this commutativity the eigenvalue functions A(w) of T(u) possess the 
same analytic properties as the Boltzmann weights (i.e. analyticity up to poles im- 
posed by the parametrization of the weights). It is well-known that a one- dimensional 
quantum spin-1 chain is associated with each 3-state two-dimensional vertex model. 
Its Hamiltonian H and momentum operator P are defined by 

tH = -^- In T| u=0 , P = -ilnT(0), (14) 

where r is an arbitrary positive scale factor. It turns out that because of the stan- 
dard initial condition ( [L3|) the Hamiltonians are sums of local terms which are given 
explicitly in the following sections. 

The plan of this paper is as follows. In section 2 we solve the models I and II. We 
calculate the partition function and the correlation lengths by employing analyticity 
of the eigenvalues A(u) and an important inversion identity. The models will turn out 
to be non-critical, independently of the spectral parameter u as all T(u) commute for 
different u. In subsection 2.3 we derive the associated quantum spin-1 Hamiltonian 
and determine its spectrum from the spectrum of T(u). 

In Section 3 the models III and IV are analyzed. They are mapped to the sym- 
metric six-vertex model and a free fermion model where the interaction parameter 
A plays the role of a crossing parameter or the chemical potential, respectively. De- 
pending on A the models show different behaviour: there are non-critical regimes 
with ferro- and antiferro-magnetic order, and a critical antiferro-magnetic regime. [] 
We also determine the spectrum of the associated one-dimensional Hamiltonians. 

Section 4 contains a discussion. 

2 Solution of models I and II 

To solve models I and II we shall apply the analytic method developed in [ 2] 
which uses functional equations and avoids the more cumbersome Bethe ansatz. We 
also refer to [2] for more detailed proofs and derivations of corresponding results. A 

*In the usual terminology of the six-vertex model one uses "ferro- (antiferro-) electric order" 
rather than "ferro- (antiferro-) magnetic order". 
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well known consequence of the Yang-Baxter equation for i?-matrices satisfying the 
standard initial condition (|13|) is the unitarity property 

J2 R h( U ) RJ i(~ U ) = VWavStf, (15) 
7,5 

where 

(p(u) = a(u)a(—u). (16) 

Explicitly we obtain 

model I: <p(u) = ~(4 - f/ 4 )(4 - i) 

model II: 9 (u) = -a* ^ - - . (17) 

Furthermore we have the additional crossing symmetry 

J2jJJ(«) = JC£{!(A-T*) (18) 

where the crossing parameter A is given in fl|). 

As shown in [2] the local relations (|i~5|) and (|18|) imply the global inversion relation 
for the transfer matrices 

T(u)T(u + A) = ipiuf [l N + 0(e~ N )] , (19) 

where 1^ is the identity matrix and 0(exp (— JV)) is a correction which is exponen- 
tially small in the thermodynamic limit (JV — > oo). (The correction term is identically 
zero for u = since the transfer matrices T(0) and T(A) reduce to right and left shift 
operators.) 

From the local crossing symmetry (|1|) and @ we also obtain 

T + (u) =T(\-u*). (20) 
Relations (|T9|) and ( p0|) directly imply functional equations for the eigenvalues A(it) 

A(u)A(u + A) = ^(u) jv [l + O^)] , (21) 
A*(«) = A(A — «*), 

which will be solved in the next subsections for the largest and next-largest eigen- 
values subject to some obvious analytical properties (and N even). We note that 
models III and IV also satisfy fllBD , but relation fll8|) does not hold. Therefore the 
method of solution for these models is different (see section 3). 
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2.1 Largest eigenvalue and partition function 

Here the eigenvalue Aq(u) which is the largest in the physical region is determined. 
It is convenient to define 



l/N 

for which we have 



*(«) = lim Ao (u) (22) 

N— >oo 



• (i) analyticity in the physical region < Rew < A, no zeros therein, 

• (ii) periodicity or antiperiodicity under u — > u + ^iri, 

• (iii) inversion identity: 

*(u)*(u + A) = ¥?(u). (23) 

For deriving the unique solution we closely follow [2] . 
Model I: The ansatz 

*(«) = \nu)F (24) 



satisfies the crossing symmetry ty(u) = \I/(A — u) for real Inserting (gj) into fl23|) 
yields a functional equation which is satisfied if only 

F(U)F(V2U) = 1-^1- (25) 



Solving this equation we obtain 



oo 



1 - 



F(U) = J] ui f +1 . (26) 

n=0 1 [/4 4 2n+2 

The function defined by fl24]) and (|26"D satisfies (i) and (ii) (it is indeed periodic). 



Therefore it is identical to the partition function per site (|22|). As explained in [2] 
the solution is unique. 

Model II: The ansatz 



leads to 



F(U)F(V^U) = 1 - — , (2* 



a 



which is solved by 

OO 1 1 



n=0 1 U A a An 



(29) 



t/ 4 a 4 "+ 2 
L • 

4n + 4 

Again the function defined by ( f27[ ) and ( p9|) satisfies (i) and is periodic in the sense 
of(ii). 

These are the final results for the partition function per site. In Section 2.3 the 
ground state energy of the related quantum spin chain is calculated from ^(u). 
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2.2 Next-largest eigenvalues and correlation length 

We now consider all eigenvalues of the transfer matrix for which 

l(u) := lim (30) 

is finite. The properties of l(u) are: 

• (i) analyticity in the physical region < Re u < A, zeros are allowed, 

• (ii) periodicity (antiperiodicity) under u — > u+^ni for even (odd) magnetization 
of the considered state, 

• (iii) inversion relation 

l(u)l{u + A) = 1, (31) 

• (iv) crossing symmetry 1*{ U ) — ^(A — u*). 

The properties (iii) and (iv) follow from (|2T|) . The magnetization for a state 
(ai, ckjv) of vertical spins (a« = 0, ±1) is defined as M = X^ =1 &i and is an even 
or odd integer. This quantity is conserved by the transfer matrix T. Property (ii) 
follows by inspecting the weights (^j), ([7]) in the product row of T (g). 

Applying (iii) twice we obtain Z(w + 2A) = l(u). Therefore, l(u) has two periods 
2A and A doubly periodic, meromorphic function is an elliptic function and is 
determined by the location of its zeros Qj and poles 0j + A. We obtain 



l(u) = Y\_ V^snh 
i=i 



AK , 

——(it 0j) 

7T 



(32) 



where snh is the elliptic function of modulus k G (0, 1) (see for instance [ 3, 4]) 
which is defined by requiring that the corresponding periods K(k),K'(k) satisfy 

K 7T 

As snh(...) is antiperiodic under u — > u + ^ni, it follows from (ii) that v is even (odd) 
for even (odd) magnetization of the corresponding eigenstate. The zeros Qj are free 
parameters restricted only by 

A 
2 

which follows from (iv). 

From the band of next-largest eigenvalues iv = 1) it is possible to derive [ 5, 
6] the correlation length as 

f = -o- (35) 



Re Qj = — , (34) 



From (|33D we find numerically £j = 308.93145 and £// = 7105.70704, for models I 
and II, respectively, which are surprisingly large values. 
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2.3 Quantum spin chains 

According to flliD there are quantum spin-1 chains associated with models I and 
II. After some calculation (see [2]) we obtain 

N 

# = E%;+i> ( 36 ) 

i=i 

with local interactions: 
Model I: 

H hJ+l = -A) - IB) - -^=(AjBj + BjAj) + ^((5|) 2 + (S| +1 ) 2 ) + 21, (37) 
Model II: 



4 12 « D 2 4 



H jj+1 = 25,- -A\-§B<- — (AjBj + BjAj) 

+ (9 - 4«)((^) 2 + (S| +1 ) 2 ) - (8 - 8a)I, (38) 

where we have used the definition 

A, = SjSj +l + SjSj + i, Bj = SjSj +1 , (39) 

in terms of the standard spin operators S x,y,z of spin one, / is the unit operator, and 
we have adjusted a scale factor r = 8/3 for model I. 

For these Hamiltonians the ground state energy per site e$ = lira E /N is easily 
calculated from ty(u) (|22|) 

reo = -(ln*)'(0), (40) 

and numerically we have (eo)j = —0.338201, (eo)n = —1.923384. 

More interesting are the low-lying energy-momentum excitations (see [2,4]) 

E-E = --(ln/) , (0) = X;e(p i ), 

P-Po = -ilnZ(0) = J> (41) 

j'=i 

with energy-momentum dispersion 

AK 



e(p) = V(l -A;) 2 + 4fcsin>. (42) 



T7T 

Obviously, the Hamiltonian has a gap A = e(0) (^ = 1) 

4K 

A = — (1 - Jfe > (43) 

with numerical values A/ = 0.0110033 and A// = 0.0018375 for models I and II, 
respectively. The decay of correlations in the groundstate of the Hamiltonian is 
described by the correlation length 



3 Solution of models III and IV 



The crossing symmetry fll8|) does not hold in the case of models III and IV. 
Therefore, the solution in these cases cannot be obtained by the methods of the 
previous chapter. However, it is straightforward to see that models III and IV can 
be mapped to the six-vertex model, i.e. a two-state vertex model with six allowed 
local spin configurations. 

The general idea of the mapping is to ignore the signs of all non-zero spin values 
in the allowed configurations of the three-state model. This is a unique prescription 
as b = and a = c for models III and IV, see (|J), (10) and (12). The remaining 
non-zero vertex configurations and corresponding weights are shown in Fig. 1. Other 
configurations than those shown in Fig. 1 are not possible because g = for models 
III and IV. As a spin 1 line could equally well be a spin —1 line, to such a line 
corresponds an "internal" degeneracy 2. 
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1 







1 1 
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a d e e p p 

Fig. 1 The allowed vertices and the Boltzmann weights of the six-vertex model of 
interacting spins a = 0,1. Note that a and d may be different (model IV). 

It is clear that the partition functions of the three-state and the related two-state 
model are identical in the thermodynamic limit. For finite systems the equivalence 
can be made correct by introducing a "seam" of modified weights in the Nth column 
of the six-vertex model 

K% - <P>K% (44) 

where (p is a "twist" which practically takes any real value. 

The equivalence of the three-state models with periodic boundary conditions and 
the six-vertex models with seam are to be understood in terms of the spectra of the 
corresponding transfer matrices acting on finite chains of length N. The action of 
the transfer matrices can be thought of as a rearrangement of zero-spins in a back- 
ground of non-zero spins without changing the sequence of non-zero spins. [Note 
that there are no intersections of lines of non-zero spins among the allowed vertex 
configurations.] Each time a non-zero spin is moved from column N to column 1 
the background configuration suffers a shift by one lattice constant which amounts 
to multiplying the initial configuration by exp (i<p), where (p is the background mo- 
mentum. In the reduced description by the equivalent six-vertex model there is no 
such change of the background configuration, but the factor exp (iip) is imposed by 
the modified boundary condition ([4*4]) . If we denote by N±,N the numbers of ±1 
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and spins within a row, which are conserved by the action of the transfer matrices 
of models III and IV, we have 



N + N + + N„ = N, 

N + — iV_ = magnetization M, (45) 

(p = multiple of- 



'N+ + N- 



In the following subsections we treat models III and IV separately as they cor- 
respond to qualitatively different symmetric and asymmetric six-vertex models, re- 
spectively. 



3.1 Model III 



Model III is defined by the requirement a = c = d, and b = g = 0, otherwise any 
values for a, p and e are allowed. The importance of the parameter A ([|) is known 
from the 6- vertex model [3, 7], and the different physical regions have different 
parametrizations (10). The associated quantum spin-1 Hamiltonian (|T4l) is easily 
calculated. It is again of the form (RBI) with 



H 



- [Aj + A 3 B 3 + BjAj) + A [(Sff + (S* +1 ) 



2B< 



(46) 



in the whole range of A with the abbreviations (|39|) . In deriving (|46|) from (10) we 
have adjusted different scale factors r in the three different regions of parametrization. 

For A > 1 the model possesses simple (frozen in) ferromagnetic order, the inter- 
vals — 1 < A < 1 and A < — 1 correspond to the critical and non-critical antiferro- 
magnetic regimes, respectively. 

For A < — 1, where A = — cosh A according to ( |10b[ ), correlation functions decay 
with length [5,6] 

f = -ss- (47) 

where k is the elliptic modulus defined by the requirement 



K' 
~K 



2A 

71 



(4* 



For — 1 < A < 1, where A = —cos A, the 6-vertex model is critical and the 
correlation functions decay algebraically. For critical systems the energy levels of 
the associated Hamiltonian ( |46|) on a finite chain of length N and periodic boundary 
conditions are expected to scale like 



TCV 



(49) 
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for the groundstate energy and like 

£ - £ = 

P = (50) 

for the excited states where v is the velocity of the elementary excitations. Systems 
which are conformally invariant [ 8, 9] have a unique groundstate and the finite- 
size amplitudes c, x, and s are identical to the central charge of the underlying field 
theory, the scaling dimensions and spins of primary fields [ 10, 11]. For instance 
the asymptotics of two-point functions are given by 

C r - ^ (51) 

From the correspondence of model III with a six-vertex model with seam [ 12, 
13] we indeed find ([49]) and (0) to hold with c = 1 and 

x = —n 2 + — — — [m - — +1 + L 

2 2(1 - A/tt) V 2W 

s = n(m- — \+ N^- + I-I. (52) 

The groundstate lies in the sector with Nq = N + + iV_ = N/2 and twist (p — 0. 
For the excited states n = No — N/2, m may take any integer value and /, J are 
non-negative integers. 

However, the numbers x and s in general do not correspond to scaling dimensions 
and spins. The problem can be traced back to the degeneracy of the groundstate 
of the system. There are about ~ 2 N / 2 background spin configurations for which 
N + + N_ = N/2 with momentum zero. (This implies a residual entropy of |ln2.) 
So, most of the correlation functions are ill-defined if they probe the background 
spin configuration, for instance (S^S*) or simply (S z ), as any magnetization between 
—N/2 and +N/2 may be realized by one of the ~ 2 N I 2 groundstates. Other corre- 
lation functions are independent of the particular groundstate like {(S^) 2 (Sj) 2 ) and 
are given by the corresponding correlation function of the six-vertex model (without 
twist). 

We summarize that only for (p — the set of exponents of algebraically decaying 
correlation functions is given by flSDD. This is a little unfortunate as for continuous de- 
pendence of x on ip we could have expected logarithmic behaviour of the correlations 
C r ~ Jdcp/r 2x ^ ~ l/(r 2x ^\nr). 

3.2 Model IV 

In this case we are dealing with an asymmetric six-vertex model as a ^ d (see 
Fig.l). The associated quantum spin-1 Hamiltonian is again of the form (EJ3) and 
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uniformly given by 



H j>j+1 = - [Aj + A J B J + BjAj] +ao[l- (SJ) 2 - (S| +1 ) 2 ] , (53) 

where the coupling a takes all real values, — oo < «o < °o- The range 1 < |«o| 
follows from the parametrization ( |12a| ) with «o = ±cothA, while the range |«o| < 1 



follows from ( |12b|) with a = ±tanhA, and we have adjusted the scale factors r = 
sinh A and cosh A in these two cases, respectively. 

To solve model IV we do not use the parametrization (12), but a different one 
which is more appropriate to the asymmetric six-vertex model (and does not satisfy 
(|13|)). The asymmetric model can be thought of as a symmetric model in an external 
magnetic field B acting on horizontal and vertical bond spins, with new weights R 
such that 

R£° = e B(°+0+v+'')ify« i (54) 

where (a/cf) 1//4 = e B . 

To satisfy fllTf ) with e = 1 the R can be parametrized by 

R\\ = R° ° = pcosv, RH = B%1 = p, R%{ — R{q — psinv (55) 

with fixed scale p = e~ 2B and an appropriate parameter v which satisfies < Re v < 
| in the physical region. The eigenvalues of the transfer matrix of the six-vertex 
model in external fields (plus "twisted boundary conditions") are given by 

where 

q(v) = Y[ sin(t> — vj) (57) 

3=1 

is a function of Nq Bethe ansatz numbers Vj. Equation ([56]) is easily derived (see 
for instance [3, 14]) as the vertical magnetic field contributes the overall factor 
q 2B(n-n ) an( j horizontal field and the twist amount to the prefactor e 2BN e lip of 
the first term on the right hand side of fl5"6]). 
Inserting (|57|) simplifies d56|) : 

A(t;) = \e 2BN+ ^ cos N v + (-l) N ° sin N v\ ft f-e~ 2B C ° S ^ ~ Vj ] ) . (58) 



The factoring is in fact due to the "free fermion" condition (|TTD . As A must be 
analytic at the zeros Vj of q(v), the first factor must vanish which determines the Vj 



t a n N Vj = (-l) No+1 e 2BN+l *, (59) 
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which means that the Bethe ansatz equations have decoupled. Inserting (|59| ) into 
fl58f) we obtain the final solution in terms of the original weights a and p 



A 

where 



1 f=i ai a ~ P®j) 



1 — p )Qj + ap 



©i = e^, fc, = ^--^. (61) 

The Jj are arbitrary, however distinct integers or half-integers for odd or even number 
iVo, respectively. Obviously, the influence of the twist angle (p is small, vanishing for 
N -> oo. 

The physical properties are most easily discussed for the associated Hamiltonian 
(|53|). Its eigenvalues follow from ( |i"4"|) 

No 

E = 2 [ao — cos fcj] — Na , 

3=1 
N 

P = J2 k J + ^ (62) 

3=1 

for both cases of parametrization (12). 

The properties depend on the coupling a . If a > 1 the groundstate is the 
iVo-vacuum, i.e. the state with iVo = 0. With a view to (|53|) this means that all 
S z = 1 or — 1 along the chain, the degeneracy is therefore 2^. For ao < — 1 the 
groundstate is the vacuum of non-zero spins. All S z along the chain are zero, i.e. 
Nq = N and all fcj-modes in ( E^j ) are occupied. The state is unique. In both cases 
we have i?o = — A^|a | and the gap to the excitations is A = 2(|aroj — 1), N — > oo. 

For | Qfo j < 1 all fc^-modes of negative energy in ( p^ ) are occupied in the ground 
state. Its energy is 

qos — N 

E = -2—^ sin -Ar - (JV - 2iVo)a . (63) 
sin jj: N 

The model is critical and the underlying conformal field theory has central charge 
c = 1. The "scaling dimensions" and "spins" (see also subsection 3.1) are given by 

(An) 2 ( ip\ 2 T - 



4 V 27T/ 

s = An(m-—)+(N-N )— + I-L (64) 

where An is the change in the particle number as compared to the ground state, the 
integer m describes an asymmetry in the momentum distribution kj, and /, I are 
non-negative integers arising from particle-hole excitations. 
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4 Conclusion 



We have studied four exactly solvable three-state vertex models and corresponding 
integrable spin-1 chains. 

For models I and II the groundstate energy of their associated Hamiltonians, the 
elementary excitations and the correlation lengths could be calculated. These systems 
turned out to be non-critical. Unfortunately, the methods of section 2 did not allow 
for the calculation of any order parameter or the multiplicity of the groundstate. 
However, the patterns of the elementary excitations indicate that model I is an 
integrable realization of a Haldane system [ 15, 16] and model II has Neel order. 
On chains with an even number N of sites the number v of elementary excitations 
may be even or odd, v — 1, 2, 3, ... . This excludes dimerization which would imply 
v even. On chains with an odd number of sites a similar analysis to that of section 2 
shows that the groundstate of model I is still separated by a gap from the first band 
of one-particle excitations, whereas model II has a groundstate which corresponds to 
the lowest edge of a one-particle band. If not conclusive, these findings are at least 
consistent with a Haldane system [ 17, 18] for model I and Neel order for model 
II. 

Models III and IV comprise different physical phases depending on the value 
of their coupling parameters A and a , respectively. The physical properties were 
derived by a mapping to the six-vertex model. We have seen that these three- 
state models contain spurious degrees of freedom which essentially do not participate 
in the dynamics. This is the reason for the high degeneracy of the ground states 
for models I and II (with residual entropy). It also is the reason why a simple 
application of conformal invariance relating finite-size corrections of energy levels to 
scaling dimensions fails. 
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